H 

r ^ 

© 

© 

fi 

n 

CO 

© 

0 

u 

r 

§ 

k i 

f 


I 

I 


DEFENSE  RESEARCH 
LABORATORY 


iiggasw 

THE  UNIVERSITY  OF  TEXAS  • AUSTIN  12.  TEXAS 


1EUTION  STATEMENT  A 

Ai-pr.  d foi  public  release; 
Dic.ribution  Unlimited 


r MOST  Project 


(/%) 

v— DRL-A-177 


ON  THE  DETERMINATION  OF  THE  RADIATION ( 
PATTERN  OF  AN  EMITTER  FROM  NEARFIELD 


by 


G.  S.y Innis^  Jr.  and  C.  W./korton 

Copy  No . : 


r. 


DRL^Acoustical  Repo^t^ 

Resulting  from  research  done  under  ^lYrrnii  nf  Hhiy 
Research  and  Development  ContractyNCJbsr-72627//' 
NE  051247-6;  NE  0511+56-4 
Problem  22-C 

Sponsored  jointly  by  the  Bureau  of  Ships 
and  the  Office  of  Naval  Research,  Code  1+11 


'D  DC 

.faEDOiffijf: 

NOV  22  1976 


7^ 


00 


u 


,a 


1"  :Q,  UTIQN  STATEMENT^ 

Appro -cd  for  public  release;  j 
Distribution  Unlimited 


DEFENSE  RESEARCH  LABORATORY 

THE  UNIVERSITY  OF  TEXAS 

AUSTIN  12.  TEXAS 


POST  GPPICK  BOX  SO20 
TCIKTYPX  AX  7S 


BOO  EAST  24TH  STREET 
GREENWOOD  S S0I2 


ERRATA  FOR  DRL  ACOUSTICAL  REPORT  NO.  177 


Page  lo,  Equation  (3.lM  The  factors  ( |2+tj2-1)  (I-tj2)^' 2 should  he 

replaced  by  (^-r,2)1/2 


Page  21,  Equation  (h.3) 
Page  22,  Equation  (5.1) 


Page  31 


Page  32 


The  term  (|  -l)  in  the  denominator  of  the 
right  member  should  read  ( g -1). 

The  factors  (l  - + i^2  - l)  should 
be  replaced  by  unity  and  m should  be  replaced 
by  zero. 

The  expansions  in  Eq.  (A.13)  are  really  for 
eikR/LjtR  rather  than  for  eiiiR/R  as  stated. 
The  statement  in  line  2 is  wrong.  One  has 
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ABSTRACT 


In  this  report  the  problem  of  determining  the  field  due  to  a distributed 
source  is  discussed.  Only  fields  in  which  the  disturbance  is  propagated  in 
accordance  with  the  classical,  linear,  scalar,  homogeneous  equation  of  wave 
motions  are  considered.  The  Green's  function  is  used  to  derive  two  solution 
forms  to  the  problem,  and  it  is  noted  that  these  same  forms  have  been  discovered 
by  other  investigators  using  different  techniques. 

One  of  the  forms  presented  requires  the  specification  of  inhomogeneous 
Cauchy  boundary  conditions  -{y—-and  ~di(r/drr4  on  the  bounding  surface  while 
the  other  form  requires  only  inhomogeneous  Dirichlet  conditions  4^4"  on  the 
boundary.  In  the  method  requiring  Cauchy  conditions,  an  approximate  relation 
between  i r~  and — dyf Oir  is  found  so  that,  in  effect,  Dirichlet  conditions 
are  specified.  The  error  introduced  by  making  this  approximation  is  mentioned. 

Finally,  the  radiation  pattern  of  an  underwater  acoustic  radiator  is 
computed  using  each  of  the  two  methods  and  is  compared  with  the  measured 
radiation  pattern. 
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CHAPTER  I 


Introduction 

The  problem  of  computing  the  values  of  a function  ! in  a closed  volume 
V in  which  ¥ satisfies  the  scalar  wave  equation  from  values  of  ¥ and/or 
c#/on  over  the  bounding  surface  S of  V has  been  studied  extensively"L,2,^,i+ 
and  there  are  apparently  three  different  methods  of  solving  the  problem. 

They  are 

1.  the  Green's  function  approach  in  which  the  solution  'S'(p)  at 
the  point  P is  expressed  as  an  integral  over  the  surface  S 
of  the  normal  derivative  of  the  appropriate  Green's  function 
multiplied  by  the  inhomogeneous  Dirichlet  boundary  conditions. 
Alternately,  the  appropriate  Green's  function  may  be  multiplied 
by  the  inhomogeneous  Neumann  boundary  conditions  and  integrated 
over  the  surface  S to  compute  the  solution  T(p).  The  use 


^Stuart  Ballantine,  "An  Operational  Proof  of  the  Wave -Potential  Theorem, 
with  Applications  to  Electromagnetic  and  Acoustic  Systems,"  The  Journal 
of  the  Franklin  Institute,  Lancaster  Press,  Inc.,  Lancaster,  Pennsylvania, 
April  1536,  vol.  221,  no.  4,  pp.  469-14-84. 

2Phi  lip  M.  Korse  and  Herman  Feshbach,  Methods  of  Theoretical  Physics, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  1953*  PP*  834-857. 


*Jaroslav  Pachner,  "On  the  Dependence  of  Directivity  Patterns  on  the  Distance 
from  the  Emitter,"  Journal  of  the  Acoustical  Society  of  America,  Lancaster, 
Pennsylvania,  January  1956,  vol.  28,  no.  1,  pp.  86-90. 


*Bevin  3.  Baker  and  E.  T.  Copson,  The  Mathematical  Theory  of  Huygens' 
Principle,  Oxford  at  the  Clarerdon  Press,  London,  England,  1950,  pp.  23- 


27. 
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of  inhomogeneous  Dirichlet  boundary  conditions  is  preferred  to  the 
Neumann  boundary  conditions  in  many  practical  situations  because 
a knowledge  of  the  function  ¥(s)  on  the  surface  S is  more  common 
than  a knowledge  of  the  normal  gradient  on  the  surface  S. 

2.  the  Kirchhoff -Helmholtz  method  in  which  the  solution  Y (p)  at 
the  point  P is  expressed  as  an  integral  over  the  surface  S of 
a function  involving  both  the  values  of  ?(S)  and  of 
°Y(S)/dn  on  the  surface  S.  The  obvious  disadvantage  here  is 
that  one  is  required  to  specify  inhomogeneous  Cauchy  boundary 
conditions . 

3.  the  modified  Pachner  method,  so  called  because  the  basic  idea 

5 

was  first  published  by  J.  Pachner.  This  method  requires 
the  specification  of  inhomogeneous  Dirichlet  conditions  over 
the  surface  S and  expresses  the  solution  ¥(P)  at  the  point  P 
in  terms  of  a doubly  infinite  series . 

It  is  the  purpose  of  this  paper  to  show  that  these  three  forms  are 
not  three  different  and  independent  solutions  to  the  problem  but  that  the 
second  and  third  forms  may  be  derived  from  the  first.  By  using  the  Green's 
function  for  a point  source  in  free  space,  one  can  derive  the  Kirchhoff- 
Helmholtz  formula  in  which  specification  of  inhomogeneous  Cauchy  boundary 
conditions  is  necessary.  By  choosing  the  Green's  function  which  vanishes 
over  the  surface  S,  one  can  derive  the  modified  Pachner  method. 


■'Pachner,  op.  cit. 


p.  86. 
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To  illustrate  the  methods  derived,  the  radiation  pattern  for  an  under- 
water acoustic  radiator  will  be  computed  from  measurements  made  in  the  near- 
field. As  the  radiator  to  be  considered  is  a long,  thin, circular  cylinder, 
the  prolate  spheroidal  coordinate  system  will  be  used.  This  choice  is  made 
because  a prolate  spheroid  can  be  made  to  fit  a portion  of  a right  circular 
cylinder  as  closely  as  one  desires.  Cylindrical  coordinates  are  avoided 
because  their  use  would  necessitate  the  specification  of  the  boundary 
values  ¥(s)  and/or  ci¥(s)/ch  over  a surface  of  infinite  extent.  One  would 
prefer  to  consider  a surface  S which  is  closed  in  a finite  region  of  the 
space  considered.  One  could,  of  course,  choose  the  surface  of  the  circular 
cylinder  contained  between  two  planes  perpendicular  to  the  axis  of  the 
cylinder;  but  then  the  surface  of  integration  (the  portion  of  the  cylinder 
plus  the  two  end  caps)  is  not  a coordinate  surface,  and  the  specification 
of  the  Green's  function  becomes  exceedingly  difficult.  Hence,  pressure 
measurements  over  the  surface  of  a prolate  spheroid  in  the  nearfield  of  the 
acoustic  radiator  will  be  used  to  compute  the  pressure  in  the  farfield 
(Fraunhoffer  zone)  of  this  radiator  using  both  the  Kirchhoff -Helmholtz 
solution  and  the  modified  Pachner  solution.  These  computed  values  will  then 
be  compared  with  observed  values. 


CHAPTER  II 


iturbances  are  propagated  in  many  media 


tude  acoustic  a 


in  accordance  with  the  classical  scalar  equation  of  wave  motions,  i.e 


containing  only  the  single  harmonic  angular 


are  of  particular  interest.  By 


equency  o>  such  that  ¥ = Ve 


,e  discussion  to  disturbance 


Equation  (2.2)  is  the  well  known  scalar  Helmholtz  equation,  and  V + k‘ 


It  can  be  shown  that  the  solution  of  eq.  (2.2)  can  be  written 


Morse  and  Feshbach,  op.  cit.,  pp 
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Expression  (2.5)  has  the  value  of  ty(r)  or  zero  according  as  r lies 
interior  or  exterior  to  the  bounding  surface  S [see  Fig.  ( 2 . 1 ) ] . The 
superscript  s denotes  values  taken  on  the  surface  S,  and  the  subscript  o 
denotes  source  point  coordinates.  n is  the  unit  outward  normal  to  the 
surf  ace  S . 

The  necessary  Green's  function  in  (2.3)  is  simply7  that  for  the 
point  source  in  free  space 


G(I^) 


in  which 


(2.4) 


s 1 

r - r 

— -o 


Using  eq..  (2.4)  in  (2.3),  one  gets  for  the  region  interior  to  the 


surface  S 


*(£>  - 5?  } L 


n F ikR  cty(r  ) 


8 / e 


R on 


ikR  1 1 


- It  — J ^o  • 


Let  only  those  source  distributions  be  considered  for  which  there 
exists  a surface  a,  which  completely  encloses  the  sources  of  the  field 
and  such  that,  for  all  points  p which  lie  on  the  surface  a,  |p|  < Q, 
where  Q is  some  positive  number.  One  would  like  an  expression  for  the 
solution  ^(0  in  the  region  exterior  to  the  closed  surface  a . To 
determine  a solution  in  this  form,  consider  the  doubly  connected  surface 


7Ibid. , p.  810. 
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a + Z such  that  all  the  sources  of  the  field  lie  in  the  volume  V bounded 
externally  by  a . Let  Z be  a sphere  of  radius  p and  be  centered  on 

Q 

the  origin.  It  is  easily  shown”"  that  the  contribution  to  eg..  (2.5)  from 
the  integration  over  Z as  p -»  °°  is  zero  provided  Sommerfieid's 
radiation  condition  ( "Ausstrahiungsbeaingung" ) holds,  i.e. 


r > , 1 

o" 

Prr  - ikv 
! op 


(2.6) 


uir* ii  o rrr. ^.y  r*6 spc c "t.  "t.o  pc^.car*  sjTi^Z-q 3 • 

Equation  (2.5)  can  thus  be  considered  as  giving  the  value  of  the 
function  y ( r ) in  the  region  of  space  exterior  to  a given  closed  surface 
S (a  in  the  preceding  paragraph)  in  terms  of  the  values  of  t(r°)  and  the 
normal  gradient  PV(r^)/6n  ; however,  n (the  outward  normal  to  a + Z) 

is  the  inward  normal  to  the  surface  S.  If  we  wish  n to  denote  the  outward 
normal  to  S,  we  write 


V(r)  - uj  j I v(S) 


ikR\ 

e 

ikR 

e 

^(s)l 

R 1 

R 

on  J 

(2.7) 


in  which  y(S)  is  written  for  y(r^)  for  simplicity.  Equation  (2.7)  is 

5 

the  well  known  Kirchhoff -Helmholtz  equation. 

The  use  of  eq.  (2.7)  for  the  calculation  of  \|r(r)  is,  in  many  cases, 
impractical  because  it  requires  a knowledge  of  Cauchy  boundary  conditions 


■"Baker  and  Copson,  op.  cit . , p.  25  (see  footnotes). 

ST. 

Ibid.,  p.  2o. 


■1 


>'hich  may  not  be  known.  In  order  to  specify  only  the  simpler  Dirichlet 
soundary  conditions,  an  expression  that  relates  v(s)  to  6y  ( S ) / on 

Ls  sought . 

For  simple  harmonic  disturbances  Y , a vector  u can  be 


w'i  Oil  y 


- " lap  sraa  v 


(2.6) 


density  p 


y is  the  acoustic  pressure  in  a medium  of  average 
is  the  particle  velocity.  The  normal  component  of  the 


t, i c -*.s  velocity  is  in  rn&gpii'fciicLs 


icop 


cy 

*5n 


(2.9) 


Unfortunately,  measurement  of  the  normal  component  of  the  velocity  is  not 


easy.  However,  if  the  direction  of  the  particle  motion  is  normal  to  the 


SUI  SlCG  o » oifl0n 
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'"'Claus  Muller,  "Electromagnetic  Radiation  Patterns  and  Sources," 

IRE  Transactions  on  Antennas  and  Propagation,  Professional  Group  on 
Antennas  and  Propagation,  July  1956,  vol. AP-U,  no.  3,  p.  229- 
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A:  ; .lev  Icr.  of  K 1 r chhof f -He Imho It 2 Formula 

The  form  of  the  solution  of  the  Helmholtz  eq.  (2.2)  in  the  prolate 
spheroidal  coordinate  system  has  been  known  since  1880.  However,  due  to 
the  complexity  of  this  form,  it  is  unfamiliar  to  many.  Therefore  Appendix 
A has  been  included  in  which  the  form  of  this  solution  and  several  of  its 
properties  are  discussed  in  detail.  The  results  obtained  in  Appendix  A 
will  be  assumed  in  the  remainder  of  the  discussion. 

The  problem  of  notation  in  this  coordinate  system  demands  a word  of 
discussion.  As  yet  no  notation  has  been  agreed  upon  by  the  various  workers 
in  the  field;  for  that  matter,  even  the  method  of  normalizing  some  of  the 

11 

functions  is  still  a metter  of  discussion.  Flamner , however,  lists 
and  discusses  the  various  notations  so  the  details  of  these  systems  will 
not  be  mentioned  here  other  than  as  pertains  to  the  notation  of  Flammer 
•which  is  used  exclusively.  Flammer' s notation  does  not  include  a simple 
symbol  for  the  expansion  coefficient  for  the  radial  functions,  but  the 
obvious  choice  for  this  symbol  has  been  made  [see  eq.  (A. 9)  of  Appendix  A] 


and  is  used  throughout  this  paper. 

Or.e  would  like  to  write  eq.  (2.13)  in  a form  such  that,  having 
measured  the  values  of  the  acoustic  pressure  \jr(s)  at  discrete  points 
on  the  surface  S,  one  can  compute  the  acoustic  pressure  at  the  field  point 


a. 


Carson  Fla 
Stanford, 


mmer,  Spheroidal  Wave  Functions 
California,  1957  > PP-  1^-15 • 


> 


Stanford  University  Press, 
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.o  tms  er.u  consider 


L R J 


; - o + ik  ' 
L R j 


ikR  ^ 
e oR 


How  if  the  point  r is  a great  distance  from  the  nearest  point  of 
the  surface  S,  i.e.,  R » 1 , the  first  term  on  the  right  of  eq.  (3.1) 
can  be  neglected  compared  with  the  second  term.  On  consideration  of  Fig  . 
(3.l)  the  factor  oR/br.  is  easily  seen  to  be 


ur.e  linds  thus 


= - cos  9 


a_  a 

on  l H J 


ik  cos  6 


for  points  r a great  distance  from  the  surface  of  integration.  Now  eq. 
(2.13)  can  be  written 


» (?)  = 


X7  /(1  + cos 


\|r(S)  dA 


At  this  point  there  are  at  least  two  paths  which  lead  to  a tractable 
numerical  approximation  to  eq.  (3.k).  Both  of  these  methods  involve  deter- 
mining an  expression  for  the  cosine  8 in  terms  of  the  field  point  and 
a general  point  on  the  surface  .As  deriving  this  expression 

involves  some  detailed,  though  not  difficult,  geometric  manipulation,  the 
first  part  of  Appendix  B is  devoted  to  this  calculation.  The  choice  of 


ORIGIN 


FIG.  3. 


I - A PORTION  OF  THE  SURFACE  S SHOWING  THE  OUTWARD 
NORMAL  n,  THE  FIELD  POINT  r,  THE  SOURCE  POINT  r 
THEIR  DIFFERENCE  R,  AND  THE  ANGLE  9 
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approaches  is  due  to  the  fact  that  there  are  probably  several  ways  of 
computing  (exp  ihR)/R  . Trfo  of  these  methods,  the  use  of  an  infinite 
series  expansion  in  terms  of  the  spheroidal  wave  functions  and  a conceptually 
simpler  geometrical  calculation,  will  be  expounded  separately. 

First,  the  more  sophisticated  method  shall  be  considered,  i.e.,  an 
infinite  series  expansion  of  (exp  ikR)/R  shall  be  used  in  eq.  (3-^) 
to  compute  the  pressure  at  the  field  point  r . With  eqs.  (A.12)  and 
(3.5)  one  can  write  immediately 


Hr) 


00  oo 

k2  V ^ 2-6om  (3) 

~2  L L ICTTT  SJ>C’^  Rmn(c'^ 


(1) 

R 

rnn 


(c,ls) 


£s  (i  - n2)*  (i  - n2)1  *1 

J L1  + - ^ + — cos*°J 


r 


r ™oU2  - D2 


(C  -O2 


s„,  (c>\)  cos  (~(v  " 90)]  v(v  5s'  ‘V  ^ 


in  which  the  surface  c = £ is  the  surface  of  integration. 

o s 

There  are  several  observations  that  can  be  made  to  simplify  eq. 

(3.5) . First,  since  the  pressure  ty(r)  at  a farfield  point  _r  is  to  be 

computed,  an  asymptotic  expansion  for  (c,|)  can  be  used  [see  eq. 

mn 


Consequently,  the  integration  with  respect  to  c?o  in  eq.  (3.5)  can  be 
performed  analytically,  and  the  field  study  can  be  restricted  to  the  plane 
9 = without  ^oss  of  generality.  Thirdly,  during  these  tests,  only 

relative  amplitudes  will  be  compared  and,  consequently,  all  constant  factors 
can  se  ignored.  These  terms  would,  of  course,  be  important  in  the  determi- 
nation of  the  absolute  pressure  in  the  farfield  from  nearfieid  measurements. 
Equation  (3-5)  then  becomes 


CO 


s=-.i  <c'r‘o)  * <V5S)  + 1o  - (1  • 


<o>* 


(3.6) 


and  to  compute  ty(_r)  at  the  point  £ all  that  is  necessary  is  that  one 
measure  v(ri0>§s)  over  the  surface  gQ  = and  apply  eq.  (3 .6).  Since 

there  is  no  dependence  on  c?o  , measurements  need  only  be  made  along  the 
portion  of  the  ellipse  determined  by  £ = and  <pQ  = C . 

As  eq.  (3 .6)  is  obviously  one  that  would  be  difficult  to  handle 
numerically,  one  considers  the  other  possibility,  i.e.,  the  determination 
of  (exp  ikR)/R  geometrically.  The  simpler  of  these  two  results  will 
then  be  used  in  this  application. 


[■ 


■ 


— . . — > - 
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As  mentioned  earlier  the  discussion  is  restricted  to  relative  values 
of  sotn  measured  and  computed  quantities]  therefore,  constant  factors  can  he 
ignored.  Also,  because  the  nearfield  data  do  not  depend  on  <p  one  can 


i'1-  - - ' 


•m  the  integration  with  respect 


to  gel 


V(r) 


.2  2 


e **  |_(l+G)  JQ  (F)  + iHJ,  (P)j  Us+VlKl-V5"  ^o’U  d\  <5' lk > 


in  which  J (F)  and  J. (F)  are  the  ordinary  Bessel  functions  of  the  first 
kind  of  orders  0 and  1 , respectively. 

In  spite  of  the  algebraic  complexity  of  the  terms  of  eq.  (5.1*0, 
defined  by  eqs.  (5*7)  to  (5. 10),  it  is  not  difficult  to  see  that  eq.  (5.1*0 
is  far  simpler  than  the  corresponding  eq.  (5.6).  The  difficulty  of  computing 
spheroidal  wave  functions  cannot  he  overemphasized.  Even  with  modern  high 
speed  computers  many  procedural  operations  and  numerical  computations 
involved  in  using  these  functions  remain  unresolved.  Simple  expressions 
and  simply  evaluated  expressions  for  these  functions  are  as  yet  unknown. 

15 

The  limited  tables  of  expansion  coefficients  of  Stratton  et  al.  required, 
after  two  years  of  programming,  ten  hours  of  computing  time  on  the  Whirlwind  I 
computer  at  the  Massachusetts  Institute  of  Technology.  Consequently,  in  the 
remainder  of  this  application  of  the  Kirchhoff -Helmholtz  formula,  eq.  (5*1^) 
will  be  used. 


" ?I.  W.  McLachlan,  Bessel  Functions  for  Engineers,  Oxford  at  the  University 
Press,  London,  England,  l$Al,  no.  50>  P-  15S- 

‘^J . A.  Stratton,  P.  M.  Morse,  L.  J.  Chu,  J.  D.  C.  Little,  and  F.  J.  Coroato, 
Spheroidal  Wave  Functions,  John  Wiley  and  Sons,  Inc.,  New  York,  19^)6. 
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Application  of  eq.  (ji.--,)  to  a situation  in  which  v(t)  ,£  ) is 

o s 

ensure  a (or  determined  in  some  other  way)  for  discrete  values  of  wil 

f c .urse , require  the  numerical  evaluation  of  the  integral  with  respect  to 

ri , ■ A thirteen  point  Gaussian  Numerical  Quadrature  fcrrr.ula  is  used  to 

evaluate  the  integral  with  respect  to  tj  . The  thirteen  point  formula  is 

o 

chose:i  more  or  less  arbitrarily.  This  number  provides  a reasonably  good 


sampling  of  V(il  »£  ) or.  the  surface  of  integration.  The  Gaussian  formula 

G S 
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is  uses  because  it  has  been  proven  to  give  the  best  approximation  to  the 
Integra-  for  a given  number  of  abscissae.  Tnis  quadrature  formula  requires 
values  *f  the  integrand  at  points  spaced  in  accordance  with  the  zercs  of 
the  Legendre  polynomials  whose  order  is  the  number  of  abscissae  to  be 
considered  (in  this  case  13). 

For  the  problem  considered  here  (the  determination  of  the  farfield 
radiation  pa* tern  from  r.earfield  measurements)  the  integrand  is  determined 
in  part  by  the  experimentally  measured  values  of  the  pressure  ^(q  , g ) 


•o' 's 


in  the  r.earfield.  Now,  if  one  considers  the  variation  of  ty ( iq  , | ) with 

O S 


r\  , one  will  find  in  general  that  \|r(r)  ,5  ) varies  not  only  in  amplitude 


o s 


also  in  phase,  i.e.,  the  value  of  ^(ti  *£  ) at  a point  (qQ,|^)  will 


have  associated  with  it  an  amplitude  p and  a phase  angle  £3.  , whereas 


at  the  point  (riT>£-)  the  amplitude  will  be  p^  and  the  phase  angle  wi! 


-e  h f 


Tnis  combination  of  amplitude  and  phase  angle  suggest  a complex 


owan,  Arnold  N.,  Davids,  Norman  and  Levenson,  Arthur,  "Table  of  the 
•arcs  of  the  Legendre  Polynomials  of  Order  l-lo  and  the  Weight  Coefficients 
.f  Gauss'  Mechanical  Quadrature  Formula,"  Bulletin  of  the  American 
Mathematical  Society,  vol.  28,  no.  ic,  October  19^2,  pp.  735-7^3- 
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representation  of  the  quantity  ,/( r>  ,1  ) as  pe*^  . It  is  in  this  sense 

' c s 

that  Y'(r,  ,i  ) is  a complex  quantity.  Obviously.  v(ri  A ) written  as 
^ b *0  s 

* 

pe*  is  actually  a complex  representation  of  a real  quantity.  Because 
of  the  concern  for  a relative  determination  of  V(r)  , only  relative  values 
of  p (to  within  a multiplicative  constant)  and  of  pi  (to  within  an 
additive  constant)  need  be  measured. 

The  values  of  p and  p3  are  measured  at  the  points  ti  =0.0 

o ’ 


± 0.225,  - 0.454,  ± 0.643,  - 0.795,  - 0.921,  and  ± O.985  or.  the  curve 
5s  - 1.0125  and  cpo  = C . These  measured  values  are  then  used  in  eq. 
(3-14)  with  the  integral  numerically  evaluated  by  Gaussian  Numerical 
Quadrature  to  determine  the  computed  radiation  pattern  shown  in  Fig.  (3 .2). 
Tr.e  measured  farfield  pattern  of  Fig.  (3-2)  is  the  pattern  for  this 
hydrophone  measured  a great  distance  (30  feet)  from  the  hydrophone. 

Tr.e  agreement  between  computed  and  measured  values  in  Fig.  (3-2) 
is  seer,  to  be  quite  good.  Other  tests  of  this  method  indicate  that  when 
the  normal  derivative  of  the  pressure  is  given  by  ik  times  the  pressure 
[as  ir.  eq.  (2.12)],  radiation  patterns  computed  from  eq.  (3*4)  are  in 
excellent  - jreement  with  observed  radiation  patterns. 


A 


!. 


r 


CHAPTER  VI 


. or..: 


Ere:;.  Figs. (y. 2)  and  ( 5 • 1 ) one  car.  see  that  both  the  Kirchhoff -Helmholtz 
: °r  -.a  a...  the  moaifieu  metho : of  Paehr.o-r  give  very  good  agreement  with 
O'-  • : va.u-jc  .r.  the  particular  example  considered  here.  Both  methods  have 

o.  vicur.y  gcoi  and  bad  points.  The  Kirchhoff  -Helmholtz  formula  car.  be 
evaluate!  ir.  such  a way  that  one  avoids  the  determination  of  the  wave 
functions,  rut  this  formula  also  requires  that  one  specify  Cauchy  boundary 
conu.tioas  even  if  only  Ly  approximation  as  in  the  example.  Tne  author  has 
attempted  to  apply  this  formula  in  cases  in  which  replacing  Cauchy  boundary 
conditions  „y  Dirichlet  boundary  conditions  via  eq.  (2.12)  was  not  valid; 
ar. as  expected,  the  computed  ar.d  measured  radiation  patterns  die.  not  agree. 

Tne  method  of  Paehr.er  requires  only  the  specification  of  Dirichlet 
cour.iary  values  but  will  ir.  many  cases  raise  more  questions  than  it  answers. 
Tne  problem  of  computing  the  necessary  wave  functions  is  a difficult  one. 

Tne  rapidity  of  convergence  of  the  series,  in  general,  and  the  dependence  of 
this  convergence  on  the  values  of  g and  j are  unknown. 

Ir.  conclusion,  the  author's  purpose  is  accomplished.  There  are  not 
- -■■'■aral  ; asically  different  solutions  to  the  problem  of  determining  the  values 
c:  a function  which  satisfies  the  Helmholtz  equation  exterior  to  a given 
..-.•face  from  values  or.  this  surface,  but  there  are  merely  several  forms  of 
one  so- .tic/..  Both  of  the  forms  of  the  solution  have  their  good  and  their 
a-  charac  eristics.  Tne  form  to  be  used  in  any  particular  situation  is 


- 2k 


APPENDIX  A 


• ’c  - :a . Wave  Function : 


i.r.e  prolate  spheroidal  coordinate  system"'1  is  related  to  the  rectangu 


Cartesian  coordinate 


syste  . via  the  systerr,  of  equations 


d 

X = 2 

: (1  - t]2)  (12 . 1)  * 

cos  9 

(A.l.a) 

d . 

y = 2 : 

. (1  - r2)  (i2  - 1)  J2 

sin  9 

(A.l.b) 

d 

2 = 2 r‘  s 

(A.l.c) 

1 « n « 1 

* ^ 5 < 00 

0 9 « 2n 

sys 


which  d is  the  separation 
The  Helmholtz  equation, 
tern,  giving  the  two  ordinary 


of  foci  in  the  coordinat 
eq.  (2.2),  is  separable 
differential  equations 


e system,  (see  Note) 
in  this  coordinate 


ar.s 


S„r(c,T)) 

toil 


2 2 
c rj 


S 

mn 


(c,n) 


c 

(A. 2) 


R (c,6) 


1 


2 2 
c £ 


0 

(A.J) 


^Fla: 


Note : 


:er, 

All 


op.  clt. , p.6 
square  roots  will  be 


taken  as  non-negative. 
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-n  which  the  m and  Xm ^ are  separation  parameters  ana  c = gka  with  k 
the  wave  number.  The  quantities  S__^(c,q)  are  the  prolate  spheroidal  angle 
functions,  and  the  R (c,0  are  the  prolate  spheroidal  radial  functions. 
Notice  that  eqs.  (A. 2)  and  (A. 3)  are,  with  a change  of  sign,  identical.  They 
are  written  separately  only  to  emphasize  the  difference  in  the  range  of  the 
two  variables  | and  q . It  should  also  be  noted  that  eqs.  (A. 2)  and 
(A. 3)  are  of  the  form 


f , 1 

— ' (I  - z*) 
az  L dzj 


L 22 

; X - C Z 


U = 0 


(A.4) 


which  is  the  well  known  Sturm -Liouvi lie  equation.  All  of  the  theory  of 
the  Sturn-Liouville  equation  is  thus  applicable  to  eqs.  (A. 2)  and  (A. 3); 
ani,  in  particular,  the  orthogonality  of  the  solutions  u in  eq.  (A. 4) 
will  be  of  interest  later  ir.  the  development. 


P'-n"'  oTo  Qt.Vi 


vr.e  r^rola- 


heroidal  Angle  Functions 


Equation  (A. 2),  a second  order  ordinary  differential  equation,  has  two 

linearly  independent  solutions;  however,  one  of  these  solutions,  the  angle 

( 2 ) 

function  of  the  second  kind  S„  ' (c,q)  , is  not  finite  for  q = ±1 

mn 

As  these  values  are  in  the  region  of  interest  in  this  particular  problem, 
the  functions  of  the  second  kind  will  not  be  included  in  the  formulation 
of  the  general  solution  and  will  not  be  considered  hereafter.  The  other 
solution  to  eq.  (A. 2),  the  angle  function  of  the  first  kind  S^,^(c,q)  , 
car.  thus  be  denoted  by  S (c,q)  without  ambiguity. 

The  solution  to  eq.  (A. 2)  is  determined  by  noticing  that  for  c = 0, 


this  differential  equation  is  the  one  satisfied  by  the  associated  Legendre 


26 


functions  of  the  first  kind  of  integral  order  and  degree.  Tnis  suggests  for 


angle  function 


ite  series  of  the  forn 


00  I 

Srr.n(c'ri)  = X dr“(c)  Cr^  (A'5) 

1 - Uj  -L 

which  the  prime  over  the  summation  sign  indicates  summation  over  only 
even  values  of  r when  n-m  is  even  and  over  only  odd  values  of  r when 
r.-m  is  odd.  The  expansion  coefficients  d'.'  ( c ) are  determined  by  substi- 
tution of  eq.  (A. 5)  into  eg..  (A. 2).  Flammer  ° gives  a recursion  relation 
for  the  d“'“(c)  and  also  finds  an  expansion  for  the  eigenvalues  X^(c ) . 

Fortunately,  however,  the  necessary  expansion  coefficients  have  been 

^ 17 

tabulated  oy  Stratton,  Morse,  Chu,  Little  and  Corbato.  1 The  limitation  of 
these  cables  to  values  of  c < 8 is  responsible  for  the  choice  of  parameters 
(specifically  frequency  and  length  of  radiator)  used  in  the  example  considered 
in  this  thesis.  With  the  tables  of  Stratton  et  al.  and  tables  of  associated 
Legendre  functions,  one  can  easily  compute  the  angle  functions  of  the  first 
kind  from  eq.  (A. 5)- 

Orthogonality  of  the  Functions  S (c,t]) 

— ■■  — ——————  mn 

From  the  general  theory  of  Sturm- Liouville  differential  equations 
it  follows  that  the  functions  S (c,rj)  are  orthogonal  in  the  interval 
(-1,1)  . The re fore  , 

/ Smn(c'ri)  Smn’  (e^}  dr‘  = &nn'  Nmn  (A'6) 

-1 


7 

‘Stratton  et 


op.  cit. 


which  the  is  easily  determined  by  using  the  normalization  factor 

.'or  the  associated  Legendre  function  to  be 


= 2 


rr.n. 


(r  + 2m)l  [d‘:n(c)] 
(2r  + 2m  + 


r=0,l 


1)  r! 


(A. 7) 


The  Prolate  Spheroidal  Radial  Functions 

As  in  the  case  of  the  angle  functions,  there  are  two  solutions  to 
the  differential  equation  satisfied  by  the  radial  functions;  and  one  of  these 
solutions  is  singular  on  the  coordinate  surface  £ = 1 . Since,  however, 

this  singularity  does  not  occur  in  the  region  of  interest  for  this  particular 
pro',  lem,  both  radial  functions  [ the  radial  functions  of  the  first  kind 

ar.d  the  radial  function  of  the  second  kind  R^^(c,|)] 

will  ■ e included  in  the  general  solution.  It  is  expedient  to  consider  the 
radial  function  of  the  third  kind  given  by 


sl5)  <c,s)  ■ ( = ,{>  ♦ (c,s) 


(A.3) 


Flammer  expands  the  function  R^^(c,|)  in  terms  of  the  expansion 

IT*' 

coefficients  d^“‘(c ) whereas  Stratton  et  _al.  introduces  a new  expansion 
coefficient  which  in  the  notation  of  Flammer  would  be  a“‘n(c)  and  which  can, 
of  course,  be  written  in  terms  of  d‘““(c).  Thus 


m oo 


R(5)  (C,|) 


mn 


L .2  J 


I 


~n  (1) 

a‘:n  (c)  h (eft) 
r ' ' m+r 


(A. 9) 


30 


(I) 

which  (cj)  is  the  spherical  Hankel  function  of  the  first  kind. 

rr.n. 


ic  expansion  coefficients  a‘"‘(c)  are  chosen  so  that 


R '(c,f;) 


_n  ^ i[c|-^(n+l)rt] 


c|- 


(A. 10) 


Expansion  (A. 9)  is  an  asymptotic  series  being  valid  only  when  the 

( 3 ) 

significant  values  of  R_  (c,|)  lie  in  the  range  0 < r < c£ 

Because  of  the  value  of  c|(c|=8.l)  used  in  the  examples  worked  out 

, ( 1 \ 
here,  another  expansion  must  be  devised  for  the  imaginary  part  of  Rw'(c,|). 

rr.n 

This  is  necessary  because  the  significant  terms  in  the  expansion  of 
( 2 \ 

(c>5)  d0  not  occur  in  the  range  0 < r < c|  . Flammer  derives  the 
expression 


i Q„n  R^(c,|)  log  + gmn  (c,5)  (A. 11) 


and  gives  details  for  computing  the  QWvi  and  g^ 


R),~ ' (e,j)  is  calculated  in  this  case  by  replacing  the  spherical  Hankel 
function  in  eq.  (A. 9)  by  the  spherical  Bessel  function.  The  radial  function 
of  the  third  kind  R^b(c>0  is  then  computed  from  eq.  (A. 3). 

Greer. ' s Function  in  Prolate  Spheroidal  Coordinates 

In  Chapter  II  the  simple  Green's  function  (exp  ikR)/R  in  prolate 
spheroidal  coordinates  is  required.  Flammer  x writes  this  as 


Flammer,  on.  clt . , p.  35- 
~^IV id . , p.  35,  eq.  (h.k.6). 

2^I_id.,  p.  35,  eq- 

" Ibid. , p.  ■;?. 


V \ 2 

2ik  / / s- 


m=0  r.=m 


_ (c,q)  S (c,n  ) cos  [rr.(cp-cp  )] 

>**»  lull  U O 


(l)  . , (5)  . , 

R_  (c,6)  R.  (c,e  ) 


5 < 5, 


(A. 12) 


(1)  (3) 

R_  (c,s  ) R),  (c,0  6>6 

■iiii  u iTi.t  O 


«’h i eh  all  terms  [such  as  N (c),  S (c,q) 


etc . j are  exactly  as  defined 


earlier.  The  subscript  o again  is  used  mo  denote  values  taken  in  the  source 

coordinates.  As  only  the  region  exterior  to  the  closed  surface  | = E 

o bs 

is  considered,  the  lower  term  of  eq.  (A. -2)  will  be  used. 

In  Chapter  IV,  however,  the  Green's  function  which  is  (exp  ikR)/R 
plus  a genera-  outgoing  wave  such  that  the  sun-,  satisfies  the  prescribed 
boundary  conditions  is  of  interest.  For  the  problem  considered  herein  the 
Green's  function  which  vanishes  over  the  surface  of  integration  will  be 
chosen  so  that  only  inhomogeneous  Dirichlet  boundary  conditions  need  be 
specified  for  the  problem.  One  gets, then, c 


G(r  rj 


+ a general  outgoing  wave 


§11 


m=G  n=m 


ir^\-  , w(5)/  Rmn  ^c,5s^  (3) 

LRmn  (c^o}  Rmr.  (c^}  “ _(3),„  t > V. 


(A. 13) 


— ' ,Js'  . v r(5)/,  , 

nrr r v c * > ^ c > * ) 


n ' o mn 


s > 6, 


22  Ibid.,  p.  eq.  (5.2.1l)-notice  that  cos  [m(<p-<p  )]  is  missing  from 
eq.  (5.2.11). 
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Ls  function  G(r|r^)  in  the  text 


|rJM0  and 


in  the  braces. 


a r w (5) 

i Rtnn 

3 Q L.  **  U Liu  I 


(D 

« (c,go)  (3) 


(c>l_) 

Kill  o 


R (c.g  ) R (c,f) 
ran  ' 3 ^o  ran  ' * 5 ' 


*3  easily  evaluated  using  the  Wronskiar.  for  (c,|)  and  R (c,|) 

(2)  (1) 

U)  ^ oR_  (c,g)  (2)  oR  (c,0 

Rrr.n  (c'^  ~ ^ \n  ~ ^ = — | • (A.14) 

5 06  cU2-l) 


(A. It)  one  gets  for  cG(rjr^)/on  [or  cG(r j r^ )/o5ol  on  the 
*o  = ^s 


°°(—  ‘-to  ) 1 


^=5s 


2rt  c(|2-l) 


I 


ra=o  n=m 


2-6 

om  „ , x 

If SmS0*') 

ran 


(3) 

, R (c,6) 

^mn  c,T^o  C0S  ^('P-'Pq)^  -7T] 

"L  (<=.6S> 


(A.15) 


Thus,  eq.  (A.13)  represents  the  Green's  function  which  satisfies  homogeneous 
Dirichlet  boundary  conditions  on  the  surface,  and  eq.  (A.15)  is  the  normal 
derivative  of  this  function. 


t 


APPENDIX  3 


To  compute  the  cosine  of  the  angle  between  the  normal  to  the  surface 


the  field  point  (5,3,0)  , consider  the  eqt 


in  the  x, z plane  (9=0  plane).  This  equation  can  be  solved 


± | (a2  - z2F 


(a  - z y 


well  known  that  dx/dz  is  the  slope  of  the  tangent  to  the  ellipse. 


sc  true  i rom  elementary  analytic  geometry  that 


m as  the  s^ope 


:e  and  m^O  , then  -l/m  is  the  slope  of  the  perpendicular  to  this 
ier.ce,  the  slope  of  the  normal  to  the  ellipse  is  given  by 


'2  . ,2^ 


a is  the  angle  the  normal  makes  with  the  z axis, 
o 


Consider  for  the  moment  the  quadrant  in  which  both  x and  z 
ire  positive.  In  this  quadrant  the  slope  of  the  normal  to  the  ellipse 


is  positive  everywhere  and  one  can  write 


1 


xpressior  f or  tne  cos  a Is  obviously  valid  for  all  x,  z 


-ipse;  novever,  tne  expression  for  tne  sin  a is  valid  onxy 

o 


per  half  plane  (x  > C)  and  must  be  replaced  by  its  negative 
wer  half  plane  (x  < 0) 

general,  the  plane  containing  the  normal  to  the  surface  of  a prolat 
and  tne  major  axis  ci  tne  ellipsoid  maxes  an  an^le  y witn  tne 
r.e . Referring  to  Fig.  (3.1)  one  sees  that  the  angle  between 


the  angle  5 : the  angle  between  r.  and  z is  0:  ; the 

o 


ween  R and  z is  a ; finally,  the  angle  between  the  plane  of 


and  the  plane  of  r.  and  z is  cpQ  One  can,  therefore, 


he  cosine  of  the  angle  9 from  the  well  known  relation  for  spherica 


cos  d = cos  a cos  a + sin  a sin  a cos  o - (3-5) 

o o o 

s merely  to  observe  that  in  the  prolate  spheroidal  coordinate  system 
major  axis  of  the  ellipse  a and  the  semi-minor  axis  b are  given 


express j y one  nns 


o U 


.u-i  )2  (1<)‘ 

w^- 


-OS  c 


'-pted  variaoo.es  denote  value: 


*e  field  point.  Notice  that  the 


ccc  vQ  tern  tarf.es  care  of  the  change  in  sign  for  sin  a mentioned 

earlier.  Thus,  the  angle  between  the  normal  at  the  point  (ri  , f , p ) 

o’  s’^o 

to  tr.e  prolate  ellipsoidal  surface  sq  = g and  the  line  joining  a farfield 
point  (q,=,C)  is  given  by  eq.  (3.8). 


To  determine  R one  can  write  simply 


a = ! (*  - *0f  + (7  - y0)  - (z  - zo)c 


2 I t 


~nd  u„e  eqs . (A.i_)  to  get  R in 


spheroidal  coordinates.  As  in  the 


Lcu-ation  of  cosine  9 , the  discussion  will  be  restricted  to  field 


r wmeh  uie  in  the  plane  <p  = 0 . Consequently,  one  gets  for 


.se  in  eq.  (3-9) 


x = | (i  - 0 (r  - Df 


y = C 


(3.10) 


z' 7 
J , 


With  the  value  of  R from  eq.  (B.lb),  (exp  ihR)/R  is  easily  found  to  be 


O O o o o U CO  O O O 


-d- 


o', 


J' 

X - 
-D-  • 


DISTRIBUTION  LIST 


C: 

Cr 

C': 


bureau  o:  Snip. 


bureau  o: 


Code  6S$A 
Code  639B 


..iei  , 


Ui. 

Of 


q v -p,  o 


Code  odO] 
Code  631 


01  Ope  1 or* s Opi»e.v  522G 


reau  of  IJaval  Weapons,  Code  AV  d3 


u: . 
Of 
Co; 


vaj  j se  arc  n . 


;ocie 


y2  IJaval  Research,  Cede  ^63 


iveiopment  Group  TV/O 


o o .... . * 1 o 


Lii: 
Di 
^ 1 
Di 

r% „ , 

v>. 


-or, 

:tor. 


rorce  Pacific  Fleet 
Officer  and  Director,  U.  S.  Havy  Electronics  Laboratory 
Officer  ar.i  Director,  U.  S.  Navy  Underwater  Sound  Lab. 
Officer,  Naval  Air  Development  Center,  Johnsville,  Pa. 

Attn : NADC  Library 

Research  Laboratory 

naooratory,  Code 


S. 


Coae  -0 


O • 


Re  sea 

c.i.  ^asoraoory, 

Research  Laboratory,  Code  2C2 
Ordnance  Laboratory 
Mission  (Navy  Staff) 

bin 

;__rer.  An  any  sis  Croup,  Brown  University 
-us.bia  University,  Lament  Geological  Observatory 
"Die  Oceanographic  Institution 
fee rat cry,  The  Johns  Hopkins  University 


--  } 

uo.:.  w fDsj 2* v . 
■?  ~ *■  q-  c • 


.recto: 


Offic 


0. 

0. 


U. 

S. 


tj  Mine  Defense  nacoratory 


u_ect.rcr.ies  Laboratory 
Attn : Dr.  Donato.  A.  Wilson 


U.  S. 


:gamo  E-ectric  Company, 
•is ion  of  Ocean: 


--  — ) 


mavy  n-ectromes  naccratory 

Attn : Dr.  G.  H.  Curl , Coue  225C 

* rvttr* : C.  H.  Lanphier 

Hydrographic  Office,  Washington  25,  D.  ( 

IP"'  '■*  q '■>  w— ^ c.  v"  ' q AMcn  ?T  V” 

-•  — - — — *-•— ev»  O*.  v*  , Oj  *.  _« C aoS  ) . a • 

Dmmittee  or.  Undersea  ’Warfare,  National  Research  Council 
rdr.ar.ee  Research  Laboratory,  The  Pennsylvania  State  University 
. C.,  U.  S.  Navy  Mine  Defense  Laboratory,  Attn:  Dr.  J.  Hagemann 
. 0.,  U.  S.  Navy  Defense  Laboratory  Attn:  Mr.  R.  J.  Urick 
oratory  of  Marine  Physics,  Yale  University 
S.  Naval  Ordnance  Test  Station,  Pasadena,  California 

Attr. : Pasadena  Annex  Library 

m-.ar.der,  Mine  Force,  U.  5.  Atlantic  Fleet,  U.  S.  Naval  Mir.ecr 

Charleston,  South 

0.  Key  West  Test  and  Evaluation  Detachment,  U.  S.  Navy  Depar 

Key  West,  Florida 

S.  Navy  Postgraduate  School,  Monterey,  California 


ui  c3  L 


'O  C\  > . \ '■  i rr\  t vo 

i \ I ill  V 


